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-  IncroTtental  angle  enclosing  material  segnent 

-  Propagation  speed  at  time  shock  has  propagated  a  distance  R. 

-  Total  kinetic  and  internal  energies  behird  the  shock  front 

-  iCinetic  energy  per  unit  mass,  internal  energy 

-  Total  nrxientum  in  material  segment 

-  i^pplied  inpulse  per  unit  curea 

-  Initial  pressure  of  gets  in  cavity 

-  Pressure  iiraediately  behind  tlie  disturbance  front 

-  Radius  of  shock  front  at  time  t 

-  Initial  cavity  radius 

-  Radius  of  cavity  during  expansion 

-  Radial  coordinate 

-  Dimensionless  cavity  radius, 

-  time 

-  Displacement  of  cavity  boundary  and  general  particle 

-  Dimeixsionlcss  shock  radius,  IVR 

-  Radial  pcurticle  velocity  behind  frornc 

-  Radial  particle  velocity  at  shock  front 

-  V/ork  done  by  the  expansion  of  perfect  gas  within  tlie  cavity'- 

-  Ratio  of  initial  to  final  specific  volume 

-  Compressibility  pairaireter,  1  -  Po/Pf 

-  Ratio  of  specific  heats  for  perfect  gas 

-  Itital  energy  deposited  on  cavity  wall,  per  unit  area 

-  Pcuranetcr  describing  energy  deposition  rate 
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An  underground  explosion  generates  a  strong  shock  wave  vAich 
propagates  tlurough  the  earth,  but  the  strength  of  this  shock  wave 
diminishes  with  distance  for  a  number  of  reasons.  In  the  case  of  a 
cylindrical  (line)  source,  or  a  sjherical  (point)  source,  gecmetric 
attenuation,  i.e.,  the  distribution  of  the  original  energy  (or  mai^- 
tum)  over  larger  material  volumes  as  the  disturioance  progresses  nec¬ 
essarily  leads  to  Icwer  shock  pressures  at  later  times.  In  this  work, 
several  models  of  tliis  process  are  proposed  and  investigated.  Waves 
diverging  in  both  cylindrical  and  sphericcil  coordinates  are  consid¬ 
ered.  Models  based  on  the  assunrotion  of  an  instantaneous  deposition 
of  a  knavn  energy  (or  mcmentum)  on  the  interior  of  a  cavity  are  devel¬ 
oped,  as  well  as  models  based  upon  a  continuous  traixsfer  of  emrgy  (or 
momentum)  from  the  material  \rLthin  the  cavity  to  the  surrounding  medium 
through  tte  inedianism  of  cavity  expansion. 

It  was  the  goal  of  this  investigation  to  develop  several  sinple 
models,  no  one  of  which  is  expected  to  be  most  appropriate  for  all 
applications.  In  any  given  situation,  the  appropriate  model  is  best 
determined  through  a  ocaiparison  of  the  various  predictions  v.’ith  that 
available  data  wliich  most  closely  notches  tlie  conditions  of  interest. 

Ti'iS.  simple  models  to  be  presented  each  satisfy  only  conservation 
of  energy  or  conservation  of  mcmentum,  rather  tlian  both,  as  the  true 
solution  must.  Consequently,  none  of  these  modeJs  should  be  regarded 
cis  e>:act  solutions  to  any  physical  problem.  In  all  cases,  tlie  material 


If 
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through  which  the  shock  propagates  is  assumed  to  be  an  ideal  locking 
iiiaterial.  Ihe  limitations  of  this  aissunption  have  been  discussed 
elsewhere  [1], 

Biis  report  is  organized  as  follows .  In  Section  II,  models  based 
on  conservation  of  momentum  are  developed.  First,  a  "snowpla^"  model 
for  an  impulse  instantaneously  applied  over  tlie  interior  of  a  spherical 
cavity  is  given.  Ihis  is  follwed  by  a  model  for  the  gradual  transfer 
of  mcii^tum  to  the  surrounding  material  by  the  expansion  of  a  spherical 
cavity  initially  filled  with  gas  at  high  pressure.  Finally,  both  of 
these  cases  are  considered  for  the  geometry  of  a  cylindrical  v/ave. 

In  Section  III  several  models  based  on  conservation  of  energy  are  dis¬ 
cussed.  First,  a  model  for  the  rapid  deposition  of  energy  over  the 
interior  of  a  spherical  cavity  is  given.  Next,  a  model  for  the  gradual 
transfer  of  energy  fran  tlie  contents  of  a  spherical  cavity  to  the 
surrounding  medium  is  developed.  Tnis  is  follo^^sd  by  tlie  analogous 
results  for  these  two  cases  in  the  cylindrical  geometry,  and  finally 
models  based  on  a  simplifying  assumption  about  the  spatial  energy  dis¬ 
tribution  are  considered.  The  report  terminates  witli  a  comparison  of  the 
predictions  of  various  models  for  a  specific  case  and  a  sliort  summary. 
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II.  ^^Clltei'lcUi'u  I''jcxlc^ls 


While  it  is  reccgnized  that  tlie  total  ncREiitum  in  a  syimietrically 
diver-^ing  cylindrical  or  sfiierical  wave  is  zero,  in  th''  absence  of 
shecir  stresses,  nvansntum  within  any  infinitesimal  solid  angle  should 
also  be  conserved.  Ghe  folla-zing  models  are  developed  in  that  spirit. 

An  extension  of  the  "sncwplOT"  model  (as  discussed  in  Ref  [1])  to 
splierical  coordinates  is  first  considered. 

a.  Prescribed  Inpulse,  Spherical  Coordinates 

We  wish  to  determine  the  rate  of  decay  of  the  pressure  pulse 
which  results  when  a  prescribed  impulse  per  unit  area  is  applied  to 
the  interior  surface  of  a  spherical  cavity  of  radius  I^.  We  assume 
that  the  material  containing  the  cavity  is  initially  of  density  p^, 
but  that  the  application  of  a  pressure  pulse  corpacts  the  material 
to  a  final  density,  p£,  v;hich  is  independent  of  the  applied  pressure. 

We  consider  a  region  of  the  material,  as  shorn  in  Figure  la,  where 
the  radius  R  is  the  location  of  the  disturbance  front  at  seme  time  after 
the  application  of  the  pressure  pulse.  and  R  are  valves  of  the 
lagrangicin  radius  coordinate  identifying  the  region  of  interest.  It 
is  convenient  to  enploy  the  radius  R  as  the  tiiiE-liJ:e  variable  in  re¬ 
cording  the  progress  of  the  event.  During  the  time  required  for  the 
disturbance  to  travel  the  distance  to  R,  the  inner  radius  has  moved  a 
distance  v/hich  is  a  function  of  tinB,  or  alternatively,  a  function 
of  the  current  disturbance  radius,  R.  At  amy  time  tlie  location  of  the 
inner  bovindary,  in  tlie  LUlerian  description,  is  as  sha-m  in  Figure  3b,  or 
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a.  REGION  OF  INTEREST 


b,  KINEMATICS  OF  CAVITY  EXPANSION 


Figure  1  Gsanetry  of  Defonnation 
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R{R)  =  Rq  +  UqCR) 

We  assume  the  splierically  syniiietric  velocity  field  has  only  a  radial 

conponent,  V^.,  and  denote  ty  V{R)  the  particle  velocity  iitriBdiately 

behind  the  disturbance  front.  In  terms  of  an  Eulerian  radius,  the 

radial  velocity  field  must  be  divergence  free  if  the  continuity  equation 

is  to  be  satisfied,  as  tlie  assumption  of  an  ideal  locking  solid  leads 

to  a  ujiiform  density  behind  the  disturbance  front.  Thus 

V„  =  V(R)  ~2  (2) 

r 

The  displaceitent  at  the  inner  cavity  radius  is 

where  t(R)  denotes  tlie  time  required  for  the  disturbance  to  reach  R. 
Denotirxg  ty  D{R)  the  rate  of  propagation  of  the  disturbance,  i.e., 

D(R)  =  ^ 

we  may  varite  Equation  3  as 


dR 

KoDTrT  [R(r)]^ 


(4) 


Substitution  into  Equation  1  yields  an  integral  equation  for  R  as  a 
function  of  R  v/hich  may  he  solved  if  V(R)  and  D(R)  are  known. 

At  the  distiurbance  front,  the  Rahldne-Hugoniot  jump  relationships 
must  be  satisfied  in  order  to  insure  tlie  conservation  of  itass  and  mo¬ 
mentum  across  the  disturbance.  Thus 


PqD(R)  =  Pf[D(R)  -  V(R)3 

(5) 

P(R)  =  pJ)[R)  V(R) 

(6) 
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where  P(R)  is  the  pressure  irrmediately  behind  the  front  measured  witii 
respect  to  the  initial  pressure.  Fran  the  first  of  these,  it  is  seen 
that  tiie  ratio  of  velocities  in  Equation  5  renains  constant  for  an 
ideal  locking  material,  and 


D{R) 

Equation  1  then  beccxnes 


=  1  -  p^Ipj  =  3 


r2 


(8) 


R(R)  =  +  6  /  „ 

R^  R(R)2 


dR 


vhicli  has  ttie  solution 


R(R).  =:[eR^  +  (1  -  3)Ro^] 


1/3 


(9) 


(10) 


This  relationship  may  be  more  easily  obtained  by  equating  the 
mass  contained  in  a  holler-/  sphere  of  density  p^  and  radii  R  and  • 
with  that  contained  in  a  holla//  sphere  of  density  p^  and  radii  R 
and  R. 


Ihe  nass  and  nrroentum  contained  in  a  length  dr  of  the  region  of 
interest  are,  fraa  Figure  lb. 
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dm  =  p£{rd(J>)  dr 

(11) 

dh  =  p^v^(r)r^dr(d<>)^ 

(12) 

At  tlae  tiiie  v^en  the  disturbance  has  readied  R, 

the  total  mcmentum 

behind  the  disturbance  is 

R 

H  =  /  o^v  (r)r''dr(d'j.)^ 

(13) 

■  R(R)'^^ 

If  an  iiipulse,  1^,  per  unit  ini.tial  area  is  sujcplied  at  t  -  0  over 

tlie  iiiner  cavity  R  ~  R^,  v/e  hcive  fran  tlie  balcuice  of  j.ripulse  and  jra-entum. 


Substituting  expression  for  the  velocity  field  (Equation  2)  into 
Equation  13,  v?e  find  that 

H  =  p^V(R)  R^  [R  -  R(R)]  {d<^)^  (15) 

where  R(R)  is  as  given  in  Equation  10. 

Canbining  Equations  7  and  8  so  as  to  eliminate  the  disturbance  speed 
yields 

P(R)  =  p  (16) 

o  e 

The  e>q>ression  for  the  attenuation  of  the  peak  pressure  as  a  function 
of  distance  is  found  by  substituting  Equations  14  and  15  into  Equation 
16,  or 

P(R)  =  ‘  (17) 

Pq  r4  [R-  R(R))^  8 

At  large  values  of  the  inner  radius  R  can  be  seen,  from 

Equation  10,  to  asynptotically  approach 

R(R)  =  R  (18) 


Hence,  at  large  R, 

P(R)  =  - -  —r/i.  2 

Po  ^  3(1  -  3^/3)  2 

i.e. ,  the  pressure  decays  as  the  sixth  paver  of  distance. 

Eor  values  of  R  only  slightly  greater  than  R^,  v/e  nay  set  R  = 
+  e,  \T^iere  e  «  I^.  Then,  from  Equation  10, 

R(R)  =  +  ae  (20) 

and  the  decay  of  pressure  with  the  distaiice  becaoes,  for  small  e. 


P(x)  =  --P 


(21) 
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When  written  in  terms  of  a  ratio  of  specific  volvcnss,  a  ,  custaiarily 


defined  through 


a  =  = 


so  Po 


=  (1  -  e) 


Equation  21  beccnes  tlie  familiar  result  of  tlie  one  din^sional,  or 

"snow-plow"  theory  [1] ,  for  the  attenuation  of  a  plane  v/avo., 

2  2 

In  a  Vso 

P  =  0 ^  (23) 

e^(a  -  1) 

vdiere  is  the  specific  volume  of  the  fully  ccrrpactsd  solid. 

Equation  17  can  be  put  into  a  form  useful  for  conparing  the 
attenuation  as  a  function  of  porosity  through  substituting  Equation 


8  into  Equation  17,  or 


P  =  P(R) 


Pf  _  (1-g)  4 

~P  8  R'^(R  - 


The  decay  of  the  dinensionless  pressure  P  is  given  as  a  function  of 
the  dimensionless  distance  in  Figure  2.  Results  for  several  values  of 
6  are  given.  For  purposes  of  ccx^>arison,  a  line  with  slope  of  -6  is 
also  sha-jn.  For  values  of  >  3  the  pressure  decays  essentially 
as  tlie  sixth  power  of  the  distance,  as  is  predicted  by  Equation  19. 

It  can  also  be  seen  that  the  predicted  rate  of  decay  appears  to  be 
slaver  for  8  =  0.9  (high  porosity)  than  for  8  =  O.i  (law  porosity) , 
but  it  should  be  noted  tliat  the  weight  of  material  required  to  produce 
a  given  pressure  reduction  will  be  substantially  less  for  tlie  material 
of  higher  porosity. 

In  Figure  3,  the  dimensionless  pressvire  is  given  as  a  function 


of  the  distance,  “  1,  measured  from  tlis  edge  of  tlie  cavity.  For 
snail  distances,  the  decay  is  seen  to  be  as  tlie  inverse  scaiaro  of  tlie 


Figure  3  Pressure  Decay  Near  Spherical  Source,  "Sna^pla-/"  I-i^del 


distance,  in  accord  witli  Equation  23,  and  as  the  inverse  sixth  po.'/er 
for  large  distances. 

Since  the  velocity  field  is  knavn  (Equation  2) ,  the  pressure 
distribution  behind  the  shodc  front  Ccin  also  be  determined.  Frcm  local 
conservation  of  nrsientuni. 


3VV  ,  3P 
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with  tile  pressure  at  r  =  R  given  by  Equation  17  arid  the  velocity  field 
given  by  Equation  2.  Ihe  time  derivative  in  Equation  25  can  be  re-' 


written  as 

_  av,  dR  _  av(R)  R-  ^  2V(R)R  J 
^t 

wiiere  the  shock  speed  D  is  also  a  function  of  the  current  sliock  radius 
R  and  the  particle  velocity  at  the  shock  front  is  given  by  Equation  15. 
Substituting  into  Equation  25  and  in-fcegrating,  we  find  that  the  pressure 


distribution  is  given  by 

P(r)  _  ^  ^  S  fl  R 

PCR)”  *  1  -  e  [8  R  -  R  \ 


)  P  1 

8  ~  11  4.  “  (L. 

R2j  r  ■  2 


In  Figures  4a,  4b,  and  4c,  the  pressure  distribution  behind  the 


front  is  given  at  several  instants  of  tiii'e  (shock  radius)  for  several 


values  of  porosity.  Ihs  substantial  grcr^th,  with  time,  of  the  inner 
radius  of  the  cavity,  K,  is  e\'ident.  Examination  of  these  pressure 
profiles  reveals  negative  pressures  in  the  region  near  the  inner  cavity 
radius.  5hese  tensions  result  frem.  the  "hoop'*  strains  becosiiing  so 
large  as  to  pro;luoe  a  net  increase  in  volume  (dilatatian)  viiicli,  for 
the  assumed  isotropic  state  of  stress,  results  in  a  radial  stress  which 


is  positive.  For  highly  porous  materials/  ha\’ever,  this  tensile  stress 
is  small  cxitpared  to  the  peak  stresses,  and  it  may  be  assumed  that  the 
solution  has  reasonable  validity  for  predicting  tlie  pressures  at  the 
shock  front  in  spite  of  the  fhysiccil  unreality  of  the  tension  near  the 
inner  cavity. 

b.  Expansion  of  a  Spherical  Cavity 

In  the  preceeding  section,  the  decay  of  the  pressure  pulse  re-~ 
suiting  from  the  instantaneous  application  of  a  prescribed  impulse 
over  the  interior  surface  of  a  spherical  cavity  was  considered.  In 
this  section,  the  decay  of  the  pressure  v?ave  generated  by  tlie  applica¬ 
tion  of  a  prescribed  pressure  will  be  treated.  This  problem  has  been 
considered  previously  and  soae  asynptotic  results  obtained  [2] . 

Ihe  expression  given  (Equation  13)  for  tlie  total  mcmentum  behind 
the  expanding  spherical  front  is  valid  independent  of  the  means  by 
which  t!ie  shock  wave  is  generated  as  is  the  kinematical  relationship 
(Equation  10)  relating  the  current  location  of  the  inner  cavity  radius 
to  the  current  shock  radius.  Hence,  if  a  pressure  P(t)  is  applied 
coniencing  at  time  zero,  the  total  inpulse  delivered  can  be  written  as 

I(t)  =/^P(t)R(t)^  (d.?)^dT  (28) 

0 

and  equated  to  the  total  rocmentum  at  that  instant,  as  given  in  Equation 
13.  ^vs  assuiis  here  that  the  pressure  supplied  is  sufficient  to  keep 
th(  .mtericil  fully  caipacted. 

Of  particular  interest  is  the  case  vhere  the  ca'-diy  is  fill©! 
vdtii  a  gas  at  high  pressure  which  expands  and  drives  a  spherical  shock 
into  the  porous  material.  Vfe  will  assume  a  perfect  gas  at  initial 
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(29) 


where  u  =  a  +  by  \  ^  v  ^  /  ■*  J 

For  Y  other  than  5/3,  Equation  34  can  be  integrated  numerically  without 

difficulty. 

5he  decay  of  pressure,  with  distance,  of  the  peak  shock  pressure 

is  given  in  Figure  5  for  y  =  5/3  and  several  values  of  3.  It  can  be 

seen  that  the  rate  of  decay  increases  witli  distance,  varying  frcra 

-3 

approximately  P/P^  -  small  to  appraximately  P/P^  - 

(R/R^)“^  for  large 

In  Figure  6,  the  decay  of  pressure  v/ith  distance  in  a  medium  of 
3  =  0.5  is  coiTpared  for  various  values  of  y  At  large  R,  the  asyitptohic 
b^avior  is  an  P  -  R  as  can  be  readily  deduced  fron  Equation  33. 

'Ihis  is  in  agreement  witli  the  results  of  Ref  [2] .  •uje  result  for 
Y  =  0  corresponds  to  a  constant  pressure. 
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Figure  5 


Decay  of  Shock\7ave  Fesultirig  fran  Expansion  of  Spherical 
Cavity,  Manentum  ^5odsl  with  y  =  5/3,  for  teterials  of  Vajrious 


Porosities 


=  0.0 


Figure  6  l>2cay  of  Splierical  Shock  Wave  in  Porous  Material,  g  =  .5, 
Kosulting  froH  Ej^nsion  of  Vcirious  Gasses,  ^fcm^ntum  I-iodel 

'7 


5.0  10.0  20.0 

R/Ro"  D5STANCE 


Y  =  1.0  corresponcJs  to  an  isotheijiial  expansion,  and  y  =  1.4  and 

Y  =  1.667  correspond  to  tlie  adia}jatic  expansion  of  air  and  a  mono  tonic 

gas,  assumed  to  be  ideal,  y  =  *5  has  no  physical  interest,  but  is 
included  for  purposes  of  coiparison.  Ihese  results  \vere  cbtained  by 
the  nunjsrical  integration  of  Equation  33  using  step  sizes  of  = 

0.001  for  1.01  £  IVHq  1.1,  =  0.01  for  1.1  £  £  2  and 

=0.1  for  2  ^  VRq  £  11.  Pressures  at  the  cjxtreme  points  of 
these  ranges  were  found  to  agree  typically  to  three  or  four  significant 
digits  whether  obtained  v/ith  10  or  100  integration  steps.  The  results 
obtained  numericeilly  for  the  case  of  6  =  0.5  and  y  =  5/3  as  given  in 
Figure  6  shew  excellent  agreement  v;ith  the  results  for  the  sane  case 
cbtained  through  evaluation  of  the  closed  form  expression  and  given  in 
Figure  5. 

c.  Mcmsntum  Models  for  Cylindrical  Waves 

Ihe  analyses  given  in  the  previous  sections  are  readily  modified 
to  be  applicable  to  an  expanding  cylindrical  wave,  as  miglit  be  generated 
by  the  detonation  of  a  line  charge  of  Mgh  explosive  imbedded  in  a 
porous  itBterial. 

In  cylindrical  coordinates  tlie  appropriate  divergence  free  velocity 
field  is 

Vr=V(R)|  (40) 

The  relationship  b8ta-;een  tlie  inner  caviiy  radius  and  tlio  shock  front 
radius  which  must  be  satisfied  in  order  to  insure  global  conservation 


or 

2  2  V2 

R  =  ler  +  (1  -  e)R^]  (42) 

which  takes  the  place  of  Equation  10,  previously  developed  for  the 
spherical  case. 

The  ncmentum  contained  in  an  element  of  angle  de,  width  dz  and 
lengtli  dr  is 

dll  =  dz  rde  dr  v^  (43) 

Thus,  v/hen  the  disturbance  has  reached  a  distance  R,  the  total  momentum 


in  the  segment  is 


R 

H  =  /  p^  dz  de  v„  rdr  =  p.  dz  de  V(R)  R(R  -  R)  (44) 

V  ^  ^  f 


If  an  impulse,  I^,  per  unit  initial  area  is  applied  over  the  inner 
cavity  radius  at  t  =  0,  the  balance  of  impulse  and  momentum  leads  to 


Pf  V(R)  R  (R  -  R)  (45) 

The  predicted  rate  of  decay  of  tlie  shock  front  pressure  then  folla^s 
from  solving  Equation  45  for  V(R)  and  substituting  the  result  into 
Equation  16 


2  2  2 
,  ^  IpC  (1  -  3) 
P(R)  =  -^-2 - 


Pq  3  I^(R  -  R)2 


(46) 


where  R  is  as  given  in  Equation  42.  This  result  may  be  termed  the 
"sna^pla-/"  model  for  the  cylindrical  wave.  At  large  R,  R(R)  -  •.  3’r 


P(R)  = 


2  2  2 

1$  R^  (1  -  3)^ 


1  1 

i  li  - 


(47) 


Indicating  that  the  pressure  decays  as  the  fourth  pcwcr  of  distance. 


at  Icurge  R. 
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As  before,  tlae  pressure  distribution  behind  tlie  sliock  may  be  found 
by  integrating  the  momantum  equation.  Bie  result  is 

fi  R  fl-B  B ')  Sfl  l£  -  l]]  (48) 

pIrT  1-e  b  ^  ^  /  r  2  \p  i 

Negative  pressures  v/ill  again  occur  near  the  cavity.  In  pcurticular,  at 

large  R,  R  6  ^'r’  and  the  pressure  on  the  surface  of  the  cavity  is 

predicted  to  be 

fflL=  1  - -i-  t  In  i  (1-g)]  (49) 

P(R)  1-8  8  2 

vhich  is  negati-ve  for  any  8. 

Ihe  solution  for  the  expansion  of  a  cylindrical  cavity  filled  v/itli 
a  perfect  gas  may  also  be  obtained  in  a  manner  similcir  to  that  given 
previously  for  the  spherical  case. 


©le  impulse  generated  by  the  expansion  is 
t 

I  =  /  P(t)  I  do  R(t)  dr  dz 
0  r=R 

Changing  the  variable  of  integration  yields 


(50) 


R 

I  =  /  P(R)  do  dz  R(0  3 


V(C) 


Equating  this  to  the  total  momentum  given  by  Equation  44  yields 


R 


fC-\  ^ 


8  /  P(R)  R(5)  vfey  =  Pf  V(R)  R  (r-R) 


(52) 


v;hicli  is  to  be  solved,  forthe  shock  front  particle  velocity,  as  a  function 
of  distance.  For  tlie  adiabatic  ejq^ansion  of  a  cylindrical  cavity'  filled 
with  a  perfect  gas. 


P(R)  =  Pq  (R^) 


2y 


(53) 


20 


Differentiating  Equation  52  with  respect  to  R,  as  before,  and  defining 


x  =  -R/R 

2 

s  =  R/R  =  tBx^  +  1-B] 
p  =  x(x~e) 

Equation  52  niay  be  written  as 


(54) 


(55) 


Replacing  by  the  pressure,  through  the  use  of  Equation  16^ leads  to 
an  expression  for  the  peak  shock  pressure  as  a  function  of  shock 
radivis 


P(x)  _  2(1-8) 


g(l-2Y)g(5_s)^g 


(56) 


Pq  X  (jts)  1 

Equation  56  bears  a  high  degree  of  reseniblance  to  the  ancilagous  expression 
previously  given  for  the  spherical  wave  (Equation  33) ,  and  was  intc.grated 
numerically  in  the  manner  previously  described.  Insults  obtained 
the  numerical  integrationof  Equation  56  are  given  in  Figure  7  for  the 
case  of  Y  =  5/3  and  various  values  of  8/  and  in  Figure  8  for  8  =  »5 
and  several  values  of  y 
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P’igvire  7  Decay  of  Cyliiidrical  Shodc  Wave  Desultj.ng  fron  Cavity 

Expansion  for  Materials  of  Various  Porosities.  ^S•:^n^nt^.Ftl 
^toc^Gl  v;ith  y  =  5/3 


III.  Energy  Models 


Ilie  models  given  in  tlie  prececding  sections  resulted  frctn  enforcing 
conservation  of  radial  iraientum.  It  is  also  possible  to  develop  predi.c- 
tions  of  tl^e  pressure  attenuation  rate  considering  the  conservation 
of  energy.  'JSie  principaJ.  difficulty  in  this  approach  is  the  uncertainty 
as  to  the  internal  energy  changes  \4rLch  can  occur  in  a  material  of 
constant  density. 

a.  Prescribed  Deposition  of  Energy  in  Spherical  Cavity 

For  a  spherically  diverging  v;ave,  the  divergence  free  velocity 
field  is 


v^ere  V(Pv)  is  the 
and  is  related  to 
bet3'/een  the  shock 


Frexn  Equations  16 


Vj.  =  V{R) 


(57) 


particle  velocity  imrediately  behind  the  shock  front 
tlie  shock  speed  by  Equation  8.  liie  kinetic  energy 
front  and  the  ej^nding  cavity  is 


R  2  r 


(58) 


cind  8,  this  may  te  rev/ritten  as 


p  -  gP(R) 
^  2(1-8) 


^-1.  (d« 


(59) 


v/here  R  is  as  given  by  Equation  10. 

Iirtnediately  behind  tlie  shock  front,  the  kinetic  energy  per  unit 
mass  is 

e^(R)  =iv(R)^ 

and,  from  the  Rankine-ilugoniot  junp  condition  for  ei^ergy,  the  increase 
in  interned  ej-iergy  per  unit  mass  is 
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(61) 


(R)  =1  /i-i 


Pf 


Substitutim  of  Equations  16  and  8  into  Equation  61  shows  that 


ei(R)  =e^(R)  =iv(R)^  (62) 

i.e. ,  that  the  internal  and  kinetic  energies  per  unit  nass  are  identical 
at  the  shock  front. 

For  a  naterial  of  the  ideal  locking  nature  under  consideration,  no 
voluna  dianges  can  occur  behind  the  shock  front.  Thus,  if  tliere  are 
no  reactions,  losses  due  to  viscosity,  or  any  heat  transfer,  the  energy 
of  eadi  particle  must  renain  constant  at  the  value  acquired  upon  passage 
through  the  shock  front.  Ihe  particle  initially  at  a  distance  r  from 
the  origin  will  therefore  acquire  and  retein  an  increase  in  energy  per 
unit  mass  of 


ej(r) 


g  P{r)| 

2Po  r=R 


(63) 


Ihe  total  internal  energy  is  then  easily  computed  in  Lagrangian 
coordinates,  and  is  found  to  be 

E  =  /  p(r)r^dr 


(64) 


If  a  Jtna.-.’n  energy  (per  unit  area)  of  is  instantaneously  deposited 
over  the  interior  of  a  spherical  cavity  of  radius  R^,  an  energy  balance 
requires  that 

e  R?(d^.)2  =  +  E_  (65) 

O  O  K  I 

throughout  tlis  entire  process.  Here  E^^  is  given  fcy  Equation  59  and 
Ej  is  given  by  Equation  64.  Satisfaction  of  this  energy  balance  for 
all  time  (i.e.,  for  all  R)  leads  to  a  decay  of  peak  shock  pressure, 
with  distance,  which  is  the  solution  of  the  equation 
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(66) 


v_  \  ^ -djaSv'VTH'ar^-^'s 


e  R^  =  ip£_fR  -  r  +  ^ 

oo  2(1-$)  !r  J  2 


■‘■f  ^  P(05^dC 


\i 


Ihis  is  a  singular  integral  equation,  the  solution  of  v^ich  is 
singular  at  R  =  for  the  physical  reason  that  a  finite  energy  is 
being  deposited  into  a  zero  volune.  In  any  plausible  physical  process, 
this  will  not  occur.  Rather,  the  deposition  of  energy  v/ill  require  a 
finite  tine.  An  arbitrary  deposition  relat5.onship 

e  =  Cq  [1  -  ei^)-  k(R-F^)/R^-]  (67) 

was  selected.  In  this  case,  90%  of  the  energy  is  deposited  in  the 
tine  required  for  the  shock  to  propagate  (£nl0)/tj  cavity  radii  into 
the  material.  Ihe  resulting  integral  equation  is 

11  -  e^-Pi  -  “W-V  A.jl  =  I  ^  [f  ' 

R  ^  2 

+  f  /  P(5)5  d5 

A  solution  to  this  equation  v;as  detained  by  nunerical  integration,  with 
results  as  given  in  Figure  9  for  several  arbitrarily  selected  values 
of  K  and  $.  Ihe  instantaneojs  deposition  of  energy  corresponds  to 

b.  E>q)ansion  of  a  Spherical  Cavity 

Instead  of  assuming  sore  arbitrary  prescripti.on  of  the  deposition 
of  energy',  we  may  assume  that  the  shock  is  generated  by  a  prescribed 
pressure  applied  over  tlae  cavity  bovindary. 

Ihe  w’ork  done  on  the  surrounding  material  during  such  an  expansion 


W=/  P(C)C^(d<;.)^dC 


where  P(C)  is  the  pressure,  a  prescribed  functicjn  of  shock  distance. 

In  particular,  if  the  pressure  results  fran  the  ej^jansion  of  a  splierical 


cavity  initially  of  radius  and  filled  v;ith  a  perfect  gais  at  initial 
pressure  P^. 

P(R)R^^  =  P  R 
o  o 

or 

w  =  [l  -  (70) 

Equating  this  work  done  to  the  kinetic  and  internal  energies  ac¬ 
quired  by  the  mterial  through  the  passing  of  the  shock  leads  to 

W  =  E^-f-Ei 

vhere  W  is  given  by  Equation  70,  E^,  1^  Equation  59  and  E^  Equation  64. 

R  I 

Hie  resulting  integral  equation  is 

(E-i,  ™ 

3(r-l)  L  R  J  2(1-8)  R  2  Ro 


This  equation  (vhich  is  not  singular)  was  solved  nuitsrically  with 
results  as  given  in  Figure  10. 
c.  Energy  Models  for  cyiindr  :.cal  VJaves 

For  a  cylindrically  diverging  v;ave,  the  divergence  free  velocity 


field  is 

Vj.  =  V(R)  I 

and  the  total  energy  behind  tlie  front  is 

2 

Ej,  ■=  ^  V^R)  ^  rdOdzdr 

^  R  2  r2 


(72) 

(73) 
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R  is  again  given  by  Equation  42.  Assuming,  as  before,  tliat  tlie  internal 
energy  of  each  particle  remains  constant  at  the  value  acquired  in  passing 
tlirougli  the  front. 


Ihus 


e(r) 


=  P(J^)6 
2P. 


R 

E  =  /  p  e(5)Cd0dzd? 


=  -  /  P(C)CdCdOdz 


(75) 


(76) 


For  a  deposition  of  total  energy,  e^,  per  unit  area,  over  tlie  interior 
of  a  cylindrical  cavity,  in  accordance  with  Equation  67,  the  energy 
balance  leads  to  a  pressure  distance  relationship  of  which  is  the 
solution  of 


2  ^ 

|l  -  exp  (-  K(R-R^)/R^j'j  =  2^1:^  P(R)R  £n(|)  +  |  /  P(C)Cdg  (77) 

Bq 

The  limiting  case  of  an  instantaneous  deposition  is  obtcdned  in  the 

limit  cis  K  A  solution  of  this  equation  obtained  by  numerical 

integration  is  sha.im  on  Figure  11,  for  several  valv-es  of  »:  and  6. 

For  the  case  of  the  expanding  cyli.nder,  the  work  done  is 

R 

W  =  /  P(R)RdRd0dz  (78) 


For  the  expansion  of  a  cylinder  of  radius  R^  of  perfect  gas  v4iich  expands 
adiabatically, 


P(R)r2y  =  p^r2y 


(79) 
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and 


W  =  |^—  'l  -  (80) 

In  order  to  satisfy  tlie  requireiiBnt  that  tlie  wrk  done  on  the  porous 
inaterial  (ICguation  80)  be  equal  to  the  change  in  kij-jetic  and  internal 
energies  (Equation  74  and  76), we  find  the  pressure  irrust  satisfy  the 


integral  equation 


R  2(1-Yr 

V  J 


&PR^  „  R ,  e 
2(l-a)  R  2 


R 

/  P(?)?df; 


(81) 


Solutions  of  this  non-singular  equation  were  obtained  nunericcd 
integration  and  are  given  on  Figure  12  for  various  values  of  B. 


d.  A  Sinple  Energy  I^odel 

The  equaliiy  of  the  internal  energy  and  the  kixietic  energy  at  the 
shock  front  noted  earlier  (Equation  62)  suggests  a  siitplifying 
approximation,  nanely,  that  the  total  energy  behind  the  front  might  be 
equally  divided  between  the  total  kinetic  energy  and  the  total  internal 
energy. 

For  the  case  of  an  instantaneous  deposition  of  energy,  it  then 
follows  that 

2  2 

^0^0^^  =  for  tlie  spherical  case 


e^R^dedz  =  2E^  for  the  cylindrical  case 

where  Ej^  is  given  by  Equations  59  and  74,  respectively. 

Ihe  results  are: 
i^’or  tlie  spherical  wave 
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(82) 
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Figure  13  Attenuation  of  Diverging  Shod^  Waves  in  Porous  Materials. 


Siaplificd  Model 
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IV.  Conparison  of  Various  Models 


A  cenparison  of  the  pressiire-distance  predictions  obtained  through 
the  use  of  several  models  is  given  in  Fig.  14  for  the  case  of  3  =  .065 
and  cylindrical  waves. 

In  each  case,  a  dimensionless,  pressure  is  presented  as  a  function 

of  distance.  In  tlie  three  cavity  expansion  irodels,  (curves  B,  D  and 

F)  tlie  pressure  is  P/Bq,  v.iiere  is  the  initial  pressure  in  a  gas 

of  Y  =  5/3.  For  tlie  "sna-;pla/'  model.  Curve  A,  the  dimensionless 
2  2 

pressure  is  the  ts-m  energy  deposition  models. 

Curves  C  and  E,  the  dimensionless  pressure  is  Pl^e^. 
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Figure  14  Precii.ctior.<5  of  Various  ^;o'Jels  for  Porous  Material  of  $  = 
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V,  Suntnary 


Several  models  for  predicting  the  rate  of  attenuation  of  the  shock 
wave  generated  by  an  ejq>losion  of  cylindrical  or  spherical  chcirges  within 
porous  naterial  have  been  proposed.  Ihe  porous  material  is,  in  all 
cases,  described  as  an  ideal  locking  solid.  Models  based  on  conservation 
of  r.>cmentuin  and  on  conservation  of  energy  are  described.  Each  model 
satisfies  conservation  of  mass  or  conserv'ation  of  momentum;  none 
satisJfy  both  and  tlierefore  all  are  to  be  regarded  as  only  approxima¬ 
tions. 

In  the  case  of  spherical  waves,  it  is  to  be  noted  that  rate  of 
decay  can  be  very  high,  ranging  from  the  inverse  sixth  paver  at  large 
distances  for  the  spherical  "snowploi';*'  nxxdel  down  to  the  third  pcs-^r 
at  large  distances  for  the  sirtplifiod  energy  models.  The  decay  rates 
for  cylindrical  v;aves  were  foujad  to  range  from  an  inverse  fourth  power 
(in  the  case  of  the  "snoi-plcw"  model)  dcwn  to  the  second  power  for  the 
simplified  energy  model. 

Fran  the  pressure  distance  relationship,  the  peaJc  particle  velocity 
can  be  determined  as  a  function  of  distance  by  using  Equation  16;  decay 
e^^onents  for  particle  velocities  always  being  one-half  tliose  cited 
above  for  pressures.  Likewise,  the  shock  velocity  can  be  determined 
as  a  function  of  distance  by  using  Equation  8,  and  can  then  be  integrated 
(nunerically)  to  give  the  shock  trajectory,  if  desired.  Ihe  time  re¬ 
quired  for  the  shock  to  propagate  to  sane  distance  R  is 


Ihe  particle  displacenent  for  a  particle  originally  a  distance  r  fron 


tlie  origin  is,  at  the  time  the  shock  hcis  propagated  a  distance  R, 


U(R,  r)  = 


|l^-(l-3)  (R^-r^)f- 


Ihe  use  of  87  would  seem  to  predict  infinite  displaosioent,  at  large 
time,  of  all  points  behind  the  shock.  Ha*;evex,  it  must  be  noted  that 
the  application  of  all  these  models  must  be  restricted  to  tiaies  (shock 
radii)  sufficiently  small  that  the  shock  speed  is  greater  than  the 
speed  of  an  acoustic  wave.  Ccsrbining  Equations  16  and  8  yields 

P(R)  =  Pq6D(R)^  (88) 
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Tlius  all  models  are  invalid  for  pressures  belw  where  is  the 

speed  of  an  acoustic  v/ave  in  the  porous  material,  and  should  be  considered 
to  be  of  questionable  validity  for  pressures  near  this  value. 
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